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1. Introduction 

The parabolic equation, originally introduced by Fock and Leontovich (ref. 1), has found successful 
application in several subject areas that deal with wave propagation issues. In the case of stochastic 
electromagnetic wave propagation through random media, the use of the parabolic equation was adopted 
by Klyatskin and Tatarskii 36 years ago (refs. 2 and 3) and has since found much success in describing 
propagation phenomena from laser beams to optical images. This type of parabolic equation is accurate 
within the paraxial approximation and is therefore suited to treat small-angle wave scattering about the 
preferential direction of propagation. The paraxial approximation holds when the wavelength A, and the 
characteristic size l 0 of the smallest inhomogeniety of the random medium is such that A«/ 0 . However, as 
shown by Klyatskin and Tatarskii, (ref. 4), this parabolic equation in the paraxial approximation is 
sufficient in describing optical propagation within the Markov approximation (i.e., where the 
inhomogeneities which compose the random medium are taken to have a 8-function correlation in the 
direction of wave propagation). There are cases, however, in which the wavelength may be on the order 
of the size of the scattering inhomogenieties A~/ 0 thus potentially violating the paraxial approximation 
and requiring an ‘extended’ parabolic equation that still possess a preferential direction of propagation but 
is not bound by the paraxial approximation, i.e., can describe scattering at large angles out of the direction 
of propagation. Such an extended equation can only be derived from the fundamental Helmholtz equation. 

It is the purpose of this work to analytically derive solutions for the generalized mutual coherence 
function (MCF), i.e., the second order moment, of a random wave field propagating through a random 
medium within the context of the extended parabolic equation. Here, ‘generalized’ connotes the 
consideration of both the transverse as well as the longitudinal second order moments (with respect to the 
direction of propagation). Such solutions will afford a comparison between the results of the parabolic 
equation within the pararaxial approximation and those of the wide-angle extended theory. In section 2, 
the extended parabolic equation for electromagnetic wave propagation is given which, of course, derives 
from the stochastic Helmholtz equation. Since this will be an operator equation in the random electric 
field, a statistical operator method is developed in section 3 which will give a general equation for an 
arbitrary spatial statistical moment of the wave field. The generality of the operator method allows one to 
obtain an expression for the second order field moment in the direction longitudinal to the direction of 
propagation. From this, expressions are obtained as a special case for the MCF of the field, both in the 
transverse and longitudinal directions. Analytical solutions to these equations are derived for the 
Kolmogorov and Tatarskii spectra of atmospheric permittivity fluctuations within the Markov 
approximation. 


NASA/TM— 2005-2 13841 


1 



2. The Wide-Angle Extended Parabolic Equation 


Consider the scalar stochastic Helmholtz equation for an electric field E(x, p) propagating principally 
along an x-axis and perpendicular to the p -plane of an otherwise arbitrary coordinate system 

^-^^^ + V 2 E{x,p)+k 2 E(x,p) = k 2 s{x,p)E{x,p) (1) 

where s(x,p) is the random part of the total permittivity s(x,p) = 1 + s(x,p) of the propagation medium. 

Using the decomposition of the total field into a forward propagating field E + (x ,p) and a backward 
propagating field E~(x, p), i.e., 


E{x,p) = E + {x,p)+ E (x,p), 


8E{x,p) _8E + {x,p) + 8E (x,p) 
dx dx dx 


with the expansion into inhomogeneous plane waves 


E ± {x,p) = | J e ± (^)exp iq -p±i[k 2 - q 2 J 12 x d 2 q 


( 2 ) 


( 3 ) 


one can obtain the following operator expressions for the component fields 

2 /^^^ ±2 (^ 2 + v ^) 1/2 i? ± (x,p) = ± k 2 (A- 2 + v 2)“ 1/2 s(x,p)[£ + (x,p)+i<-(x,p)] (4) 

such an equation (or its Fourier transform) was considered by Malakov and Saichev (ref. 5) as well as by 
Klyatskin (ref. 6, p. 169) and Frankenthal and Beran (ref. 7). 

In the event that the backscattered field is insignificant with respect to the forward scattered field, i.e., 
E + (x,p)» E~(x, p), one can formally set £' _ (x,p) = 0 ; equation (4) thus becomes a single equation for 
the forward propagating field 

2 i dE+ ( x ’^ + 2 (k 2 + Vp ) 1/2 E + (x,p) = k 2 [k 2 + V 2 )' 1/2 s (x, p )E+ (. x , p) (5) 

This equation, which can be called the “extended parabolic equation” (ref. 2, p. 169), was previously 
analyzed by Saichev (ref. 8) using a method different from the operator formulation in what is to follow. 
Its integral formulation is also known as the method of multiple forward scatter (refs. 9 and 10). However, 
in (ref. 10) it was pointed out that this method is applicable to situations where the wavelength X and the 
smallest scale of inhomogeneity l 0 is such that ).</(, (ref. 11); The classical parabolic equation in the 
paraxial approximation holds for cases where 7_<</o- Hence, although equation (5) is capable of 
describing propagation situations in which the wave is scattered at angles up to nil with respect to the 
x-axis, its application is limited to cases where ax/q. 
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The form of equation (5) can be simplified by defining the differential operator 


f V 2 
Up =2k 2 l + -± 


and the corresponding integral operator 


1 f V 2 

v =—— i + _p 

p “2 kA k 2 


Equation (5) then becomes 

2ik 8E ( X ’& + U p E(x, p) - 2k%e(x, p)E(x, p) = 0 (8) 

dx 

Hereafter, the superscript ‘+’ will be dropped. This equation can be reduced to the well-known 
stochastic parabolic equation in the paraxial approximation by expanding the operators to give 
Up ~ 2k 2 + V ( 2 and V p ~ l/ 2 k 2 and transforming the field via E(x,p)~ W(x,p)cxp(ikx) . In what is to 

follow, this relation will be employed to yield an operator equation for the generalized statistical moments 
of the field. Here, the term ‘generalized’ connotes moments at differing transverse coordinates p as well 
as differing longitudinal coordinates x. 


3. Operator Solutions for the Generalized Field Moments 

3.1 An Expression for the Generalized nm th Field Moments 


Defining the stochastic operator 


D xp = 2ik^—-Up +2Ar 4 F p s(x,p) 
ox 


equation (8) simply becomes 


D x ,pE{x, p) = 0 


which is easily amenable to further statistical analysis. To this end, one defines the generalized moment of 
the electric field 

t " 1 nm (■*• 1 jP 1>^-2 jP 2 >‘ "XmPn j-'-n+t sPn+lj" '^n+m >Pn+m ) = ( Sum ) ( 11 ) 


gnm=\\ E (xj^Pj) Y\E*{x h pi) 

j = 1 l=n+\ 
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Using a modification of the prescription set fourth in, e.g., (refs. 12 and 13), one can employ 
equation (10) for each of the field points (x,p) and obtain for the product of fields 


n§nm 




,X 2 ,P2,---X„,pn,X n+ i,Pn + l ," ' X n +m’Pn+m 


) = 0 


( 13 ) 


where 


L 


nm 


n n+m ^ 

2 ik y—+ z — 


+ 


Z (^p, + 2/f 4 Fp ; e(xy,p y ))- 


y=i 


w+m 

- Z fe + 2* 4 F p * rfa.p,)) 

l=n+\ 


(14) 


In order to isolate the quantity (g,„„) = T mn from this relation, it is expedient to adopt the methods of 
(refs. 14 and 15) and decompose the operator L nm into its average and random parts, i.e., 


( 


(/- 1 nm ) — 


Y— ■ + 

.,dxi 

j = i 


Z 

l=n+\ 


A 

dx 


7=1 


Z<4 


P j A-i Pi 

/=«+! 


(15) 


and 


Lyim — 


ri n -r rn 

Z V p,Z ( x j>Pj)~ Z ^P/£*(u,P/) 


7=1 


/=;; + ! 


= 0 


(16) 


Hence, equation (13) becomes 



Ensemble averaging this relation yields 


(j-'nm )l _ n 


\ L n m§nm } ^ 


Similarly writing 


Snm ~ f nm E Snm > (Snm ) — 0 


and substituting into equation (18) gives 



0 


(17) 


(18) 


(19) 


(20) 


Remembering that it is the goal of this development to obtain an expression for the general field 
moment V nm , one follows the development given in reference 15 and subtracts equation (20) from 
equation (13) and using equation (19) obtains 
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Combining the first and third members of this equation using the fact that 
\j-‘nm — (A?hj —LnmS-nm gives 


LnmSnm T ^nm^nm y-‘nm§nm j ~ 0 


(22) 


One must now isolate the random quantity g nm by defining an operator L n ) n inverse to L nm , i.e., 
L~jnL„ m - 1 . Thus, operating on equation (22) with L~}„ yields 


Snm ^nm^nm^nm L nm y-'nmSnm J ~ 0 


(23) 


Finally, operating on this relation with L nm , ensemble averaging and solving the resulting expression 


for ( L nm g nm ) gives 





(24) 


Substituting this result back into equation (20), one obtains for the equation governing T nm 



I 


>r 


T T~ 1 T 
^nm^nm^nm 



= o 


(25) 


The solution of this operator equation gives an exact solution for the arbitrary field moments for 
wide-angle propagation through a random medium characterized by the stochastic permittivity s(x, p) 

and the assumption that (L nm ^ = 0 . Just as in the case of the application of the local method of small 

perturbations (ref. 12), an explicit assumption governing the statistics of the fluctuations s(x,p) does not 
need to be made in the evaluation of the ensemble averages, as will now be demonstrated. 

The general relation given by equation (25) can be reduced to the parabolic equation for the field 
moments in the paraxial approximation in the case where /.«/. In this instance, one employs the 
approximations for the operators U p ~ 2k 2 + Vjj and V p ~ 1 / 2k 2 used earlier. In addition, the classical 
parabolic equation considers statistical moments in the same transverse plane, i.e., Xj = X/ = x . Thus, the 

partial differential operators in equation (15) collapse into the single operator d/dx. Equations (15) and 
(16) then become 


2 *- + 2 ^- 2 *; 


3 = l 


P j Pi 

l=l+n 


(26) 


and 
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( 27 ) 


r =lr 2 

J-'nm. — 


n n+m 

+2]s(x,p / )- 

j = 1 / = 1 + H 


Two related palatable approximations must now be made; since L nm = (L nm ) + L nm , and it is usually 


assumed that L nm « 1 , one has 


m ~ (L nm ) L nm ] ~ (Lnm ) 5 


/-I ~ 
L-'nm ~ 


1 \L n mLnm ) I ~ 1 


r 


(28) 


where 


-l 


^ n n+m 

2il-+tv 2 - Tv 2 

dx +i p J Pl 

j = 1 /=!+« 


2ik — 
dx 


-l 


2/A: 


— f dx' 

ik J 


(29) 


Equation (25) then becomes 


n+m 


2 zl ^+Z V p / - X V p, ~(L nm (L-} n )L n 

j = 1 / = ! + « 


r =o 


(30) 


where 


Lnm ( L,im )L nm ) ~ 


A " 3 

2/ 


3 .a. 

rJ 


n u u 


i= 1 £=1 


i=l /=!+« 


»+w n+/« 


- X X (s* (ac.py ) e(jc', pi )) + X X ( 8 *( x ’py) g *( x '’p 0 ' 

j=\+nk=l j=l+nl=l+n 


dx' 


(31) 


which is the well-known classical paraxial form for the problem (ref. 13). It is interesting to note that the 
‘geometrical optics’ approximation made in equation (29) leads to the parabolic equation in the paraxial 
approximation. Thus, one can envision a substantial extension of this development beyond that of the 
classical treatment if one is to use the entire form of the operator i.e., use the inverse of the 

operator (L ltm ) as solved in the paraxial approximation rather than in the geometrical optics approximation 
which was used above. This will form the subject of the next section in which the first order statistical 
moment and the generalized second moment (i.e., the mutual coherence function) are derived using the 
above formalism. 
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3.2 Solution for the First Moment 


The first order moment of the random electric field in a plane transverse to the direction of 
propagation is defined through equation (1 1) to be given by 


ri 0 (x;pi) = (£(x;pi)) 

which is a solution of the operator relation of equation (25), viz, 

|( L lo)~ l-^io^fo)]" (Ao^io^T 


(32) 


to 7 r x to 


r, n =o 


(33) 


where, from equations (15) and (16), 

(i l0 ) = 2ik — + Up ,Z 10 -2^% 1 8(x,p),p = p 1 (34) 

The ability to proceed in an analytical fashion is dependent upon the simplifying approximations used 
earlier, viz, one has that 


Ao - '(Ao) + Ao]~ ~(Ao) \ 


1 - 


( Ao ^)]" 1 


(35) 


Using this result, equation (33) becomes 

j(Ao)~^Ao(Ao) ^10^10 =o 

Employing the appropriate definitions of the operators, equation (36) gives 
2 ikj- + U p -/(2k% e(jc,p))|2iA:^ + C/p| (2£ 4 U p s(x,p))) 


(36) 


T 10 (jc,p) = 0 (37) 


This differential equation in the operators U p and V p must now be simplified and solved for the first- 
order moment T|o(x,p) . 

To this end, one must first deal with the factor 


-l 


2 ik — + U p [ = G(x, p) 

dx 


which is the Green function of the operators 2 ik d/dx + U p , defined by 

\g(x, p) = S(x - x')8(p - p') 


(38) 


3 -* 

( V C>] 

1/2' 

2 ik 1 2 k 2 

1+^- 


dx 




( 39 ) 
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where the definition of U p is used. Applying the approximation U p ~ 2 k 2 + V 2 , and solving for the green 
function G(x,p) yields. 


G(x, p) = G(x, p; x\ p') = f — ) exp[- ik{x - x ')]^ ■ ^ P \ /2 ^ x x ^ - (40) 

x — x' 


Thus, the third term within the brackets of equation (37) can be written 

-l 


(•••) = ((2£ 4 Tp &{x,p^\lik-^- + U p \ (2k 4 F p s(x,p)) 


= (2 k 4 )“J J G(x,p;x',p') (Fp's(x,p)Tp'8(x,p'))(7 2 p'(&' 


(41) 


0 —00 


Proceeding further, one now must deal with the operator products 


Fps(x,p) = 


_lYi + iT 1/2 

2 k 2 k 2 j 


:(x,p) 


(42) 


Since s(x,p) is a random function, it can be represented in the form of a Fourier-Stieltjes integral 
(ref. 16), i.e., 


;(x, p) = | exp(/k • p)/z(x, k) 


(43) 


in which the spectral amplitude dz{x, p) is endowed with the same statistical properties as is the random 
function s(x,p) as will be shown in what is to follow. Applying equation (43) to equation (42) gives 


Tps(x, p) = - J 2 exp (* ■ pV Z ( x ’ *) 


(44) 


Thus, the ensemble averaged product appearing in right side of equation (4 1 ) becomes, 

,2 


(f p 8(x,p)^s(x',p')) = ^j 

• exp(/k • p + /k' • p') (dZ{x, k )c/Z(x',k')) 


2 \~V 2 f K '2Y 1 ' 2 
1--V 

V k J 


(45) 


One now makes use of the fact that the atmospheric permittivity fluctuation field s(x,p) is taken to 
be statistically homogeneous, characterized by a power spectral density ® E (x,k) in the transverse plane, 
and 5-correlated in the longitudinal direction; these circumstances allow one to write (ref. 1 6) 


(dZ{x, K)dZ{x\ k')^> = §(k + k')F E (x-x',k)i/ 2 Ki/ 2 K' 


(46a) 
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where for S-correlated fluctuations in the x direction, the two-dimensional spectrum F z 
by 

F g (x - x', k) = 2nd(x - x')O e (k) 


(x-x',k) is given 
(46b) 


in which ®(k) is the three-dimensional spectrum of permittivity fluctuations. Using these relations in 
equation (45) and performing the integrations where possible yields 


/ 1 \2 oo , 2 \ 1 

(u p 8(x,p)Fp'8(x',p')) = 27r[^— 8(x-x') J ll-p- exp[nc-p rf 


]® e (kV 2 k 


(47) 


where p^ = p - p' is the difference coordinate. 

Taking the statistics governing the random field s(x,p) to be also isotropic, i.e.,® E (x,k) = ® e (x,k), 
equation (41) can now finally be evaluated by substituting into it equations (40) and (47); converting the 
integration in the p^ -plane into one in plane polar coordinates and performing the associated integrations 
gives 


(•••) = -/ti 2 /c 3 | 
o 



® e (k)ki/k 


(48) 


where the 8-function relation 


A 1 

J 8(x - x')dx' = — (49) 

o 2 


is employed. 

Returning to equation (37) and, substituting equation (48) into equation (37) gives 


2 ik 1-2 k 2 

dx 


r v 2 ^ 1/2 

1 + 71R 
k 2 


+ in 


K 


2h 


inV \ l‘-F 


-1 


® e (k)ki/k 


Tio(x,p) = 0 


In the plane-wave case, one has that 


f ^ V/2 


, + 5i 

k 2 


r io(*>p) = 


vf /2 

k 2 


rio(x)=ri 0 (x) 


(50) 


(51) 


since the plane wave will not possess any transverse variations. In this special case, equation (50) 
becomes 


2 ik — + 
dx 


2 k 2 +in 2 k 3 j" 



-1 

i H) 
0 v J 

® e (k)ki/k 


Tio(x) = 0 


(52) 
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the solution of which is 


rio(x) = rio(o)exp 


ikx -^—k 2 xj 
0 


'.-iT" 

. 


® e (k)ki/k 


(53) 


This result differs by the factor (1-k 2/ k 2 )" 1 from that of the parabolic equation in the paraxial 
approximation. Of course, the later is obtained from the former by retaining the first term in the series 
expansion of the factor. The presence of this factor tends to accentuate the spatial frequencies near the 
value of the wave number k. Since this theory is applicable to those situations in which /.</o, most of the 
contribution of this factor to wave scattering will occur at the largest spatial frequencies of the 
inhomogeneities. Hence, in the case of atmospheric turbulence, one can employ the von Karman spectral 
density 


O £ (k) = 0.033C|(k 2 +^ 0 2 )“ 11/6 (54) 

which is not bounded at the high spatial frequencies, to compare the result of equation (53) with that of 
the paraxial approximation. Substituting equation (54) into equation (53) and evaluating the integral 
(ref. 13) yields 


Tio(x) = Tio(o)exp 


ikx - \ (l .772 + 1 .023/)/c~ 5 / 3 + 0.3908^o 5/3 2^1 


1 K 2 

1 T — • 2- 

”6’ k 2 


\C 2 k 2 x 


(55) 


where the structure parameters for the permittivity and refractive index are related by C| = 4C 2 . In the 
case of the open atmosphere, one always has within the bounds of the extended parabolic equation, 
k»K 0 ; thus, the hypergeometric function reduces to unity and k 5I3 «K 0 ~ 513 , allowing equation (55) to be 
approximated by 


r 10 (x)=rio(o)exp 


ikx-0390SK 0 5/3 C 2 k 2 j 


(56) 


which is the result of the parabolic equation in the paraxial approximation. Hence, the use of the extended 
parabolic equation only makes negligible amplitude and phase corrections to the first order moment of the 
wave field. This result establishes the accuracy of the parabolic equation in the paraxial approximation for 
the first order moment (mean field) as it applies to atmospheric turbulence. The next section will consider 
the calculation of the generalized second-order moment (generalized mutual coherence function) of the 
wave field from the extended parabolic equation and compare its result to that of the paraxial 
approximation. In addition, due to the completeness of the operator analysis, one naturally obtains 
expressions for the MCF along the longitudinal axis. 

3.3 Solution for the Generalized MCF (second moment) 

The generalized MCF of the random electric field in two planes transverse to the direction of 
propagation is, from equation (11), given by 

rii(xi,pi;x 2 ,P2) = (£ , (xi;pi)£'*(x 2 ;p2)) (57) 
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which is a solution of equation (25), in this case given by 

[(At>-[l- (A t A - ! 1 )Y (a 1 - Ai 1 A 1 )}a 1 = 0 


( 58 ) 


where, from the definitions of equations (15) and (16), 


(Al) s 2 ik ^ + 2 ik^- + U pi - u; 2 

OX i OX 2 


(59) 


and 


A 1 - 2 k 4 [v p] e(x,,pi)- F p * 2 s* (x 2 , p 2 )J (a i) = 0 


(60) 


As with the case for the first order moment, two related approximations must be made at the outset to 
render the problem analytically tractable. In particular, so long as [A 1 1 « 1 , 


A l — [(A l ) + A l ]~ ~ (A l ) j _1 (A l A i ) 


(61) 


which allows equation (58) to become 


(Ai)-^Ai(Ai) 1 Ai^ jrii — 0 


(62) 


At this point, it is suggested to connect the longitudinal coordinates Xi and x 2 to the related centroid x c 
and difference Xd coordinates, 


Xl + x 2 

x c = — — — , X d =X\~X 2 


(63) 


the operator expressions of equations (59) and (60) then become 


(L n ) = 2ik — + U pi -U* 2 


(64) 


Ai =2k 4 


A>, 8 


Xd 


-Pi 


T/ r* 

~Vp 2 z 


x d _ 
X C - — ,P2 


Ai) = o 


(65) 


Hence, equation (62) can be written as 






2Ar 4 (H pi 8 x c H ~,Pi -V P2 s* x c -^f,p 2 


P2 




K*"'' 


-v 


2^ 4 r Pl 8 




pi 


-^2 8 


, X d _ 

*c ~~-P2 




An =o 


( 66 ) 
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The solution of this equation commences with obtaining an expression for the Green function 


2 ik— + U pi -U;A =G{x c ,p u p 2 ) 


Proceeding as in the last section and using the approximation U p ~ 2k - + V p , this requires the 
solution of 

{ 2ik l^ r + Upl -t/ p2| G! ( x ’Pi’P2) = 5 ( x c-4)8(pi-p0 8 (p2-P2) (68) 


which is given by 


( ik\( 1 1 

G(x c , pi,p 2 ) = - — — 7 - exp [/A- ((p ! - p ' 1 ) - (p 2 -p' 2 ))l2{x c -x c )] 

v ^ J\2.n) \ x c x c 


= G(x c -x ' c , pi -p'i,P 2 -P 2 ) 


Therefore, the fourth term in equation ( 66 ) is given by 


(•••) = (2£ 4 ) - j j j g{x, p\,p 2 \x',p[,p' 2 )(y Pl v p[ &* x c -—^-,p\ s x;+^-,pi 


0 -00 -00 v v ' 

- V Pl V P 2 *( X c + : y.Pi z*ix' c - ^-,P2 -V* 2 V p ,s*fx c -^-,p 2 zix' c +^-,p'l + 


T/ r* Xd — 1 X d 

+ V P2 V P2 8 ^"’P2 8 7 

V ^ J \ ^ 


,p 2 8 *^-^->P 2 ^jd 2 p'\d 2 p' 2 dx' c 


One now employs the Fourier-Stieltjes transform as before to represent the products of VpS, i.e., 


V p E^x c ±^,pi t2 l ~^ 2 L ex p(^i, 2 -Pi, 2 ^^c±^,Ki >2 ,etc. 


and obtains the following relations (functional arguments have been suppressed but the correspondence to 
those in equation (70) follows) 


/ j \2 / ^2 A 1 

V Px V p'^z) = T7T J !— 7J exp[/K, - (pi -p'\)}F z {x c -x c ,k\)d 2 K { 


r pi r p\ 88 *) = (j^r] ex p[«i -(pi -p' 2 )]F z {x c - x' c + x d ,k\)dh 
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V P 2 V P [ S * S ) = ex P[-^2-(p2-p'l)] F e(^c-4-^ ) K 2 V 2 K2 (72c) 


Fp* 2 Fp^s*s*) = f-^yl jfl-||- exp[- nc 2 ■ (p 2 - P 2 M* c - 4,k 2 V 2 k 2 


Substituting equations (69) and (72a) to (72d) into equation (70), performing the required integrals 
over p/ and P 2 ' i n plane polar coordinates, and taking all the spectra F E (x,ki >2 )= F e (x,ki >2 ), i.e., to be 
isotropic in the frequency k , one obtains 




- exp i 7, ~ -'o ( K P</ ) F e fa -x' c + x d ,K)- 
2k 

-exp -i K %- J 0 (Kp rf )F E (x c -x;-x rf ,K) + 


+ exp i 


%7 — •A>( k P«/Ve(*c -x;,k)1k£/k£/x, 


where p^ = Pi -p 2 is the difference coordinate. Using equation (46b) in equation (73) and performing 
the x c integration, remembering equation (49), finally gives 


(•••) = -2 Tt 2 /k 3 jfl-^- 1 1 - exp ^(kp^I^Wki/k 


Hence, equation (66) becomes 


2ik — + U pi -U; 2 + 


271 2 z1 3 J I-JY 1 1 - exp - 1K Xd x/ 0 (kp d ) iO E (k )ki/k I) 1 = 0 (75) 


where r n = T n (x c jCd, p</). Making the plane wave approximation of equation (75), analogous to that done 
earlier for equation (50), the resulting differential equation has as a solution 


1 1 (x c ,x d , Prf ) = r, 1 (0,0, p rf )exp - n 2 k 2 x c j 1-^- 1 1-exp ^o(kPp/)U £ (k)ki/k (76) 


(More explicitly, one has in equation (75) 
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(^ Pl -^p 2 )rii(^i,pi;^2,p2)~(v^ -v2 2 jr 11 (x 1 ,p 1 ;x 2 ,p 2 ) 

= 2V* •Vp rf r 11 (xi,x 2 ,p rf ) 

-0 

where i? = (pi + p 2 )/2 is the centroid coordinate.) It is important to note the initial condition T^O.p,/); 

since one necessarily must take x c , one must also have x d = 0 by equation (63) since ,V| ,x 2 >0. 

This result cannot be analytically studied in its entirety and hence will be considered in two special 
cases. The first case is defined by x d = 0 in which one deals with the transverse MCF 
r(x c ,p rf ) = r (x„0,pj. Thus, the extended parabolic equation solution for the MCF in a transverse plane at 
a distance x c from the source is, from equation (76) 


f 1 1 ( x c ' Pc/ ) = ft l (0,0, Pd )exp 



{l _ *7o ( K Pc/ )] < 1 ) e (k)kc/k 


(77) 


As with the case of the first-order moment, this result differs from that of the paraxial approximation 
in the presence of the factor ( 1 -icf /rf 1 which serves to accentuate spectral contributions for frequencies 
near k. The integral indicated in this expression cannot be analytically evaluated using the von Karman 
spectrum, equation (54). However, since the integrand is such that no singularities exist in the use of the 
unbounded Kolmogorov spectrum, viz., equation (54) with K 0 = 0, one can use such a spectrum in 
equation (77) and, upon evaluating the integral, obtain 


ril(* c -Pc/) = fll(0,pc/)exp 


- (3.544 - 2.048/>- 5 /3(l - J 0 {kp d )) + 


+ 0.108/t 2 p ll /3 1 /7 2 


1 ? 17. * 2 PT 

S ’ s' 9 . 

V 6 6 4 y 


\k 2 C^x c 


(78) 


where jF 2 (---) is a generalized hypergeometric function and, as noted earlier, Cf = 4C ( ) . At the outset, 

since /.«p d in most applications, the first term within the braces of equation (78), although quite 
interesting in structure, is negligible with respect to the second term. The hypergeometric function of the 
second term is most easily dealt with by first converting it to a Lommel function s^oC") 
(ref. 17, p. 986), i.e., 


I r 2 


1. 17 17. k 2 p 2 d ) 
’ 6 ’ 6 ’ 4 


= "T ( k P d) 


1 - 11/3 


%.o 


(kp d ) 


(79) 


The Lommel function reduces (ref. 17, p. 985), in the case where kp d »\, to the simple approximate 
result sy Q {kpd ) ~ (kp d ) 5,/3 . Using this in equation (79) and (78) becomes the paraxial result 


~i \ ( x c->Pd)= l”i i(0,Py)exp \A51k 2 C 2 p 5 J 3 i 


(80) 
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Again, as with the first-order moment, corrections to the second moment afforded by the extended 
parabolic equation are negligible in the case of atmospheric turbulence. This is due to the fact that the 
Kolmogorov spectral density level near the inner scale of turbulence is much smaller than it is at larger 
scale sizes. If another spectral density were considered, e.g., one describing propagation through an 
aerosol medium, this may not be the case and significant corrections to the transverse MCF may prevail. 

The second special case in which equation (76) will be examined is in the instance where the factor 
(1-K 2 /k 2 )~ 1 can be neglected; as shown above, this is a good approximation for atmospheric turbulence. 
Thus, one is now dealing with the solution for the generalized MCF for field locations at different 
transverse and longitudinal points, 


T u {xc,x d ,p d )=T u (0,0, p d )exp 


- U~K~X C J 
0 




Ik 


( 81 ) 


Longitudinal correlations of the wave field have been previously considered in (refs. 18 and 19) using 
different methods and obtaining different results. Setting x d = 0 in equation (81) gives the well known 
paraxial result for the transverse MCF T(x c ,0,p</). Unlike the cases studied above, the form of 
equation (81) suggests the use of the Kolmogorov spectrum as modified by Tatarskii, which incoiporates 
a cutoff at high spatial frequencies by allowing the introduction of the inner scale of turbulence / 0 , viz., 


® e (k) = 0.033C|k n / 3 exp 


' K 2 " 




i^m — ' 


5.92 


(82) 


Substituting equation (82) into equation (81) and evaluating the resulting integral yields for the 
generalized MCF 


Tt l (* c , x d , p d ) = Ti i (0,0, p d )exp 


- A352k 2 C^x c \B 5 ^ \F\ 


_5 P T 

6 ” 4 B 


- K 


-5/3 


(83) 


where 


B = -L + i*L 

k£ 2k 


(84) 


Thus, as first noted in (ref. 19), the presence of the diffraction factor on equation (81) modifies the 

I 1 1/2 

effect of the cutoff frequency k,„. In the case where p d > 1/k^ + ix d /2/cl ' , equation (83) reduces to 

equation (80) upon employing the asymptotic representation of the confluent hypergeometric function 
i F \(. . .). When p d = 0, equation (83) becomes 



x d j 0)= Tn(0,0,prf)exp 



f ■ 2 A 

5/6 1 


-4.352k 2 C,7X c K,„ 5/,3 < 


-1 

> 


l 2k J 




( 85 ) 


In the appropriate limits, this expression gives 
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( 86 ) 


A i (x c , x d ,0) = Ti i (0,0, p d H 


1 .8 1 'i'MkC^n x c x d _ xfk « l o 
- (0.563 1 + 2A02i)k 1 / 6 C^x c x 5 J 6 \x d X » 1$ 


Hence, there is a phase variation in the longitudinal direction, as expected, with an attendant 
attenuation as the longitudinal separation is increased. It must be noted, however, that one needs to realize 
the condition Xrf>/o in all cases so as to satisfy the assumption of 8-correlation of the fluctuations along the 
longitudinal axis (ref. 19). Figures 1 and 2 show plots of equation (85) with rn(0,0,0) = 1 for a typical 
atmospheric propagation scenario at X = 0.63 pm for two values of inner scale size, l 0 = 1.0 mm and 
lo = 1.0 cm, respectively. 
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rii 


X = .63 pirn, Inner Scale=l . 0 mm 



Xd 


Figure 1. — T 1 1 (x c ,x,j,0) versus x d (in meters) 
for A = 0.63 mm, l 0 - 1.0 mm, x c = 5 km, 
and C£ = 10“ 12 m“ 2/3 
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X. = . 63 jjim ^ Inner Scale=l . 0 cm 



Figure 2. — T n (x c jc d ,Q)) versus x d (in meters) 
for X = 0.63 mm, / 0 = 1 .0 cm, x c = 5 km, 
and C 2 = 10 -12 
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